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lO ■ I- INTRODUCTION 

o ■ 

I J/-0,77c mesons production in e+e~ annihilation at energy = 10.6 GeV remains very interesting task for theo- 

^ . retical investigations. The cross section of the process e'^e" — > ipr]c measured at Belle experiment was first presented 

in paper Q . The lower bound of the cross section measured at Belle 



X 



(T(e^e ^ ipT]c) > 33fb 



is about an order of magnitude higher, than the theoretical predictionsy] obtained in the framework of NRQCD 0. 
fSl Some efforts were made to explain this discrepancy. For example, in [j| it was assumed that some of Belle's J/tprjc 
J> [ signal could actually be J/tpJ/tp events but later in it was shown that QCD corrections decrease the value of 
'/"J ' e+e" J /tpJ /tp cross section and subsequent Belle analysis excluded this possibility completely. Among the 
I other possible explanations the contributions from glueball or color-octet states. See also a complete review of the 
Quarkonium Working Group |9| and refrences therein. 

There was a hope for the improvement of theoretical prediction by higher order QCD corrections. As was shown 
ly-j in paper QCD corrections really increase the cross section by a factor 1.8, but this is still insufficient to reach 
• experimental results obtained at Belle. 



Recently a surprisingly simple solution of this problem was found. It turns out that taking into account the intrinsic 
^ ^ motion of quarks inside charmonium mesons one can significantly increase the value of cross section. First this effect 
I was observed in the framework of light cone expansion method in works In paper this effect was considered 

as the expansion of the amplitude in relative velocity of quark inside mesons. 
^ ' In the last paper it was proved that NRQCD series for the amplitude of the process e+e~ — > iprjc in relative velocity 
• • of quark-ai itiq uark pair in the ip, rjc mesons converges slowly (in potential models relative velocity for these mesons 
IS f ~ 0.5 [IJ). The reason of such behavior consists in strong dependence of quark and gluon propagators in the 
diagrams of the process from relative momentum of quark-antiquark pair in meson. In order to show that this effect 
really takes place for the process e"'"e~ —^ iprjc the amplitude was expanded in relative velocity series except the 
I propagators of intermediate particles for which the exact expression was used. As the result NRQCD prediction is 
multiplied by the factor that regards internal motion of quark-antiquark pair inside mesons and the cross section 
becomes in 2-5 times greater than that in the framework of NRQCD depending on the width of the wave function 
used. Thus one can conclude that the usage of the leading approximation of NRQCD proved to be unreliable for 
e+e" — > iprjc at energy ^/s = lO.GGeV. 

As it was mentioned already the other method that can be used for theoretical prediction of cross section of 
the process e~ — > ^?7c is light cone expansion method. The cross section calculated in the framework of light 
cone does not contradict to the experiment data. Unfortunately in papers 0, ^| this method was used only for 
the calculation of the e~ tprjc cross section. It would be interesting to look how this method works in other 
reactions measured at the experiments. For instance in addition to the process e"'"e~ — > ipi^c Belle collaboration 
has measured the cross sections of the processes e+e" ip{2S)'nr., ipnJ^S), il){2S)rir.{2S),ipYr n. ^ f25'')Yrn[l5l|. Lately 
BaBar experiment has measured the cross sections of the processes e+e~ — > ip-q^ 'ip'Hci^S) ,%pXcoU^- In the frame work 
of NRQCD these processes were considered in paperQ. As in the case of e+e~ ^r]c Belle and BaBar results are 
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in contradiction with NRQCD predictions. In our paper we will consider leading order contribution to the processes 
e+e~ ip{2S)r]c,tpr]c{2S),^p{2S)r]c{2S) in the framework of hght cone. As to the processes e+e" VXcOi V'(2<S')xcO 
we will argue that the leading order contribution to the cross section in 1/s series is much less than NLO. We are 
going to calculate NLO contribution in our forthcoming publication. 

This paper is organized as follows. Next section is devoted to the consideration of the processes e+e^ 
ip{2S)ric,'ipT]c{2S),ip{2S)ric{2S) in the framework of light cone. In section 3 we argue that light cone LO contri- 
bution to the cross sections e+e^ — > t/'XcO, '0(2S')xcO is much less than NLO one. In last section we summarize our 
results. 



II. THE PROCESSES e+e" 'iPric,i'{2S)'nc,TpVc{2S),ip{2S)nc{2S). 

In this section the processes e+e" VP, where V — ip{lS),ip{2S), P = r]c{lS),ric{2S) will be considered. The 
cross section of the process can be written as follows 

<e-e-^ VP) = (1) 

where -^i is invariant mass of e+e~ system, p is the momentum of final mesons in the center mass frame. The 
formfactor F^^ is defined in the following way 

(y(pi,A),P(p2)|JMlO) = W-e'^Kp^i^.p. (2) 
The asymptotic behavior of formfactor F^^ with mesons in the final state can be obtained from the following formula[T^ 

(M(pi,Ai)M(p2,A2)|J^|0) ~ f-^j . (3) 

Obviously in the case M(pi,Ai) = ^(pi, A), M(p2, A2) — P{vi) the hehcity A2 = 0. As to the helicity Ai it is seen 
from formula (0) that vector meson V is transversely polarized Ai = ±1. So the asymptotic behavior of the amplitude 
is 

(F(pi,A),P(p2)|J,jO) ^ (4) 

or Fy^ ~ 1/s^ is the asymptotic behavior of the formfactor. 

Two diagrams that give contribution to the amplitude of the processes under consideration are presented in Fig 1. 
The other two diagrams can be obtained from the depicted ones by charge conjugation. The leading order contribution 
to the formfactor was first obtained in papersJTj, .12] where the process e+e~ i\)r]r. was considered. In our paper 
we follow work [l^l . In deriving the expression for the formfactor F^^ in paper the mass difference of the final 
mesons '0 and r]c. was disregarded. The mass difference of ■)/; and rye mesons is about ~ lOOMeV and this value can not 
give large correction to the cross section. But if, for instance, ip{2S) and r]c mesons are considered the mass difference 
is about ~ 700MeV. As it will be seen from the subsequent analysis this value is large enough to give large correction 
to the cross section under consideration. So we have derived the formula for the formfactor taking into the account 
different masses of final mesons. The expression for the formfactor Fyp can be written as follows 
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Where ~ (s — My — Mp), Pa, Pp,Vt,Vl,V±,Va are the hght cone wave functions defined in Tl, Mv,Mp 
are the mass of the vector and pseudoscalar mesons correspondingly, Mq = Mq^ = Mq^), Zt and Zp are the 
renormaUzation factors of the local tensor and pseudoscalar currents, d{x, y), s{x), s{y) are defined as follows: 
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Zmip^) = 



Zm ~ Z„i{k'^), Z^-^ ~ Z„j(ct^) , 



where Mq{ii^) is the running MS'-mass, Cf — 4/3, bo = 25/3, k = {ki + h) is virtual gluon momentum, a — 
—ki — li — I2 is virtual quark momentum in Fig. lb. For the light cone wave functions Pa, fp, Vr, Vl, V±,Va of 15* 
state mesons the following expressions will be used: 



4>i{x, 



c,{v')<t>t{x) 
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[1 - AXIX2{1 - v^)] 



(6) 



where w is a characteristic speed of quark- antiquark pair in meson, Ci is the coefficient which is fixed by the wave 
function normalization j dx4>i{x,v'^) = 1, is the asymptotic expression for the wave function 

In order get the wave functions of 2S states the following procedure will be used. We recall that 25" state Coulomb 
wave function has the form: 



^2s{r) (1 - qar) exp(-go7-) = 1 + go 



dqz 



(7) 



where go — qs/'^ is the mean momentum of quark inside meson, ^'is(r) is the 15 Coulomb wave function with Born 
momentum equals go- In momentum space 25* wave function has the form: 



l + go- 



1 



dqoJip^ + qlY ip^ + ql) 



- 3gg 

2\2 



(8) 



2S wave function © has zero at = iq^. Obviously this zero has nothing to do with real zero of cc{2S) meson 
wave function. In order to connect wave function (jS]) with more realistic model we replace 1 + qo-£^ by 1 + /Sqo-^- 
The constant f3 is fixed by the condition that zero of the modified wave function must coincide with zero obtained 
from the solution of Schrodinger equation with potential [14( . Thus we obtain /3 = 0.38. Now it is easy to find wave 
function of 25 state: 



i + /3go 



X1X2 



dqoj I [1- 4x1x2(1 -w^)] 



ax\X2 



[1 - 4x1X2(1 




X1X2 



[1- 4x1X2(1 



(9) 



The constant a equals unity in the case of usual Coulomb wave function. In paper the constant a was taken 
1 — since this value allows one to link different behavior of the wave function: w — > and t; — > 1. In our analysis 

we will take the same value of this constant. Finally one gets the light cone wave functions Pa, Pp, Vt,Vl,V±, Va of 
25 state mesons 



(1 — U^)xiX2 



X1X2 



[1 - 4x1x2(1 - v^)] / [1 - 4x1x2(1 - V^)] 



(10) 



The asymptotic expression for the wave function (j)'^ are given as follows 
for the leading twist 2 wave functions: 



PAix) = Vl{x) = Vrix) = r{x) = 6x1x2 , 



(11) 
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for the non-leading twist 3 wave functions: 

Pp(x) = l, V^{x)^-^[l + {x^-X2f]. (12) 

Va{x) = Pt{x) = 6x1X2 . 

The expression for Ught cone wave function ® , pU|l is one of the possible ways to link different limits: quark-antiquark 
pair in meson being in the rest — > and very light quark v ^ 1. In the former limit one obviously gets 5{x— 1/2), 
the later one leads to the function ~ 0°. 

In the numerical analysis the following parameters will be used: 

Mc = 1.2GeV, 
^(15), 7^,(15) |/p| ~ ~ 0.41GeV, 

^(25),7y,(25) |/p| ~ ~0.28GeV, (13) 
The values /y was obtained from decay width T{V — > e^e") 

r(^-e+e-) = i^^. (14) 
^ ^ 27 My ^ ' 

The constants /y, /p are considered to be equal: fp^fv- For a sip.) one loop result will be used 

47r 

oo log(M /A ) 

with A = 200MeV. L ast p arameter needed for numerical analysis is the width of wave function v'^. It will be taken 
from potential models p^: 

■iP,T]c = 0.23, 

V'(2S'),77c(2S') = 0.29. (16) 

The result of the calculation is presented in Table 1. The second and the third columns contain experimental result 
measured at Babar and Belle experiments. In the fourth column the results of this section are presented. In order 
to compare the result with NRQCD predictions for the processes under consideration the fifth column contains the 
predictions in the framework of this model. 

From table 1 one sees that the predictions of the cross section of the processes e+e" 
4>ric,ip{2S)ric,4'Vc{'^S),'4'{'2'S)i]c{2S) in the framework of light cone is much greater than NRQCD predictions. As 
was noted above the reason of this discrepancy may be attributed to the fact that at leading approximation NRQCD 
does not regard the motion inside final mesons. In paper 0] it was noted that NRQCD corrections to the amplitude 
with final tp{2S),r]c{2S) mesons are large(expansion parameter in this case is about ^ 0.7). So the apphcation of 
NRQCD to this processes is unreliable. Moreover strong dependence of the amplitude from the propagators of the 
intermediate particles mentioned above does not improve NRQCD also. As the result NRQCD prediction is in poor 
agreement with the experiment data. In contrast to NRQCD the leading order light cone predictions is in better 
agreement with data, from what one may suppose that light cone expansion is more reliable for e"*"e~ VP at 
energy ^/s = 10.6GeV. The problem with light cone expansion is connected with poor knowledge of light cone wave 
function (Pl. lfTUIl especially in the case of 25 meson. And in order to get better understanding of the processes under 
consideration in addition to QCD corrections and next to leading term in 1/s expansion one should obtain better 
knowledge of wave functions if^. ljlOII . 

There are two contributions to the formfactor Fyp separated in formula The first contribution originates from 
the wave function of mesons at the origin and it is proportional to ~ fvfp- The second contribution regards internal 
motion of quark-antiquark pair inside mesons and it is proportional to Iq. As was noted above leading NRQCD 
approximation does not take into the account contribution of the second type. Moving form the lower cc states to 
the upper ones we diminish the value of the constant fv, fp- So in the framework of NRQCD the cross sections for 
the production of upper cc mesons is less than that for the lower cc states. This effect is well seen in table 1. The 
second contribution ^ /q where the internal motion is taken into account compensate first effect since upper lying 
resonances are broader. To the first sight one may conclude that the cross section of the processes e"'"e~ il){2S)rjc 
and e+e^ — > ipr]c{2S) can not differ significantly. But our calculations show that the cross section e+e^ ipr]c{2S) 
is almost two times larger than that for e"^e~ ip{2S)ric. The reason of such large discrepancy consists in the fact 
that some terms in formulajSjis multiplied by the factor Mp/My- This factor enhances or diminishes different terms 
in (0 what results in the enhancement of the cr(e^e~ %Ij7]c{2S)) in comparison with cr(e+e^ tl>{2S)'qc). This fact 
is a peculiarity of light cone prediction. 
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III. THE PROCESSES e+e" ^ iP(lS)xco,i'{2S)xco- 

Leading asymptotic behavior of matrix element {V{pi, X), S{p2)\Jfj,\Q) may be derived from formula Q. Here we 
have M(pi,Ai) = V{pi, X), M{p2, X2) = S{p2)- Obviously the heHcity A2 equals zero. As to the vector meson V the 
leading contribution is given by the helicity Ai = 0. So the asymptotic behavior of the amplitude is 



{Vip^,X),Sip2)\J,\0) ^ (17) 



Then the asymptotic behavior of the cross section a{e^e^ VS) is ^ Unfortunately leading in 1/s expansion 

contribution is much less than NLO. To prove this NRQCD result for the cross section of the processes under 
consideration obtained in paper will be used. Let's consider the process e~^e~ — > •ip{lS)xcO- Cross section of this 
processes can be represented in the form 



^3 



a = l^a^alqlF.r^VY^^MA. (18) 

where = 16m^/s and Fq = 2(18r-^ — 7r*)^ + r-^(4 + lOr^ — ir^y. Let us substitute s — > 10.6^t and expand the 
above formula in 1/t series (toc — 1.4GeV). We get 

. = H5+yf+o(i/«»). (19) 

Thus one sees that in the framework of NRQCD NLO correction at energy -^s = 10.6 GeV is an order of magnitude 
larger than that leading one. In light cone NRQCD result is multiplied by the factor that accounts internal motion. 
Thus if one supposes that these factors are of the same order of magnitude for LO and NLO contributions one may 
conclude that NLO contribution in the framework of light cone is much larger than LO. The same is true for the 
process e^e~ ip{2S)xcO- Leading order result can be found in paper [18|. NLO contribution will be considered in 
our forthcoming publication. 



IV. DISCUSSION. 



In this paper we have reanalyzed the production of excited charmonium mesons pair in e~^e annihilation (i.e. 
the reactions e+e' 'ilj{lS)'nc{lS),'ilj{2S)T]c{lS),'ip{lS)rici2S),ip{2S)7]c{2S) and VXco) in the framework of the light 
cone. It is shown, that the internal motion of quarks inside charmonium mesons leads to substantial increase of the 
cross sections of these processes and reasonable agreement between theoretical predictions and available experimental 
data can be reached. It is also shown, that at energy = 10.6GeV NLO contribution to the cross sections e+e^ 
'iJj{lS)xcOi^{'^S)xcO a-re about an order of magnitude higher, than the LO contribution. So in order to achieve 
agreement between theoretical predictions of the cross sections of the processes e^e~ 'iP{IS)xcOt'4'{'2S)xcO and 
experimental data one should calculate NLO contribution to the cross sections. 
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